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Abstract 

We introduce a modification of the Navier-Stokes equation that has the remark- 
able property of possessing an infinite number of conserved quantities in the 
inviscid limit. This new equation is studied numerically and turbulence proper- 
ties are analyzed concerning energy spectra and scaling of structure functions. 
The dissipative structures arising in this new equation are curled vortex sheets 
contrary to vortex tubes arising in Navier-Stokes turbulence. The numerically 
calculated scaling of structure functions is compared with a phcnomenological 
model based on the She-Leveque approach. 

Finally, for this equation wc demonstrate global wcll-posedncss for suffi- 
ciently smooth initial conditions in the periodic case and in R 3 . The key feature 
is the availability of an additional estimate which shows that the L -norm of 
the velocity field remains finite. 

Keywords: PDEs, turbulence, Navier-Stokes equation, Burgers equation, 
regularity 



1. Introduction 

In this paper, we introduce a new equation which is a hybrid of the Navier- 
Stokes equation and the Burgers equation. This equation possesses a nonlocal 
part in the nonlincarity analogous to the pressure term in the Navier-Stokes 
equation. This new equation has the remarkable feature of possessing an infinite 
number of conserved quantities in the inviscid limit. 

Turbulence properties of this equation are analyzed using numerical sim- 
ulations. We calculate energy spectra and scaling of higher order structure 
functions. A key observation from the numerical simulations of this modified 
Navier-Stokes equation is that the most dissipative structures consist of curled 
vortex sheets contrary to vortex tubes in conventional Navier-Stokes turbulence. 
Using this information, a She-Leveque type model is derived and compared to 
the numerically obtained scaling of higher order structure functions. 



Email address: grauerOtpl.rub.de (Rainer Grauer) 



Preprint submitted to Physica D 



March 4, 2013 



In addition, for this new equation we can show existence and regularity for 
H 1 initial conditions of arbitrary size. This will be carried out in R 3 and periodic 
domains in R 3 . The simple modification of the nonlincarity makes the proof of 
global solutions possible, insofar as an additional estimate is available showing 
that the solution remains finite in L p , 2 < p < oo. With p = 4, this is then 
coupled with standard estimates for the -f/^-norm to complete the proof. 

The outline of this paper is as follows: in section 2 we motivate and intro- 
duce our new equation. Turbulence statistics and phenomenological modeling is 
considered in section 3. Section 4 contains the proof of existence of global solu- 
tions. We finish with remarks on possible further consequences of the existence 
of an infinite number of conserved quantities. 



2. Model equation 

We consider a three-dimensional domain f2 which shall be either R 3 or a 
bounded cube in R 3 with periodic boundary conditions. Let P = 1 — A -1 V <g> V 
be the Leray-Hopf projection operator (with periodic boundary conditions when 
is bounded): 

P[P[u]] = P[u] , V • P[u] = . (1) 
The usual incompressible Navier- Stokes equation 
d 

—v + v • Vv + Vp = vAv + f , V • v = 

at 

can be written with the projection operator P 

— v + P[v- Vv] =^Av + P[f] , V-v = 

ot 

such that no explicit pressure term is present in the equation. 

We can rewrite the Navier- Stokes equation without the incompressibility 
constraint in the form 

d 

— u + P[u] • VP[u] = i/Au + f , (2) 

where the solution of the Navier-Stokes equation can be recovered by taking 
v = P[u]. 

The equation ^ can be compared with Burgers equation whose structure 
is formally similar: 

d 

— u + u Vu = ^Au + f . (3) 

at 

For equation the nonlincarity is purely local, whereas for equation @ the 
nonlinear interaction involves the nonlocal projection. 

A natural hybrid of these two equations leads a new model equation involving 
a compressible velocity field u that is convected by its solenoidal part P [u] : 

d 

— u + P[u] • Vu = vAu + f . (4) 
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More accurately this means: The convection of the velocity field u is local in 
position space, but the projection operator is local in Fourier space and thus 
shares this mixture of local and non-local interactions with the original Navier- 
Stokcs equation. 

3. Turbulence statistics 

By construction the presented model equation is an intermediate step be- 
tween the Navier-Stokes and Burgers equation, which in turn differ significantly 
in their dynamical evolution and turbulent behavior. In Navier-Stokes turbu- 
lence, on the one hand, the most dissipative structures are vortex filaments, 
while for Burgers equation shocks dominate the turbulent flow. It is of obvious 
interest in how far our model equation bridges between those, which structures 
are the most dominant for turbulent flows and how these structures influence 
the turbulence statistics. We therefore extend the She-Leveque reasoning, which 
describes Navier-Stokes and Burgers turbulence well, to our model equation and 
test it against numerical simulations by comparing the scaling exponents of the 
structure functions. 

Numerical simulations are carried out with a second-order in space finite dif- 
ference scheme with a strongly stable third-order Runge-Kutta time integration 
with resolutions up to 512 3 . The initial conditions were chosen as Orszag- Tang- 
like (see 0]) large-scale perturbations: 

u x = A (— 2 sin(2y) + sin(z) + 2 cos(2y) + cos(z)) 
u y = A (— 2 sin(x) + sin(z) + 2 cos(x) + cos(z)) 
u z = A (sin(x) + sin(y) — 2 cos(2x) + cos(y)) . 

For simplicity and comparability both velocity and its solenoidal projection 
are set to equal values. The physical domain stretches from — tt to tt; the 
above defined conditions, thus, are both large-scale perturbations and periodic. 
All hydrodynamical models will be simulated in comparison, using these initial 
conditions. We consider only decaying turbulence without external forces. For 
the parameters of all performed runs see Table [1] which shows the numerical 
value of the quantities at the time of maximum enstrophy t = tg. 

FigureQ]sliows the decay of kinetic energy for the considered hydrodynamical 
models in comparison. The tendency of Burgers turbulence to form shocks and 
the dissipative nature of these structures lead to a faster energy decay compared 
to the Navier-Stokes equation. The new model equation exhibits a less violent 
form of dissipation, its energy decay lies in between the others. The difference 
in turbulence development is identified in a more precise way when comparing 
the time ts of maximum enstrophy £ = J n ui 2 dx. As Figured] (right) indicates, 
the enstrophy of Burgers turbulence reaches its peak significantly faster than for 
the Navier-Stokes equation, in which vortex filaments dominate the turbulent 
flow. The proposed model equation ranges between them. This hints at the 
development of coherent structures at a timescale slower than shock-formation 
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PDE 


N 


Ax 


te 


V 




£ 


L 


V 


Rx 


Burgers 


128 
256 
512 


0.0491 
0.0245 
0.0123 


0.461 
0.603 
0.634 


0.0241 
0.015 
0.009 


1.761 
1.637 
1.656 


3.564 
3.320 
3.361 


1.533 
1.320 
1.352 


0.0445 
0.0318 
0.0216 


40.99 
46.49 
61.10 


Model 
equation 


128 
256 
512 


0.0491 
0.0245 
0.0123 


1.438 
1.459 
1.624 


0.01 

0.006 
0.0036 


1.316 
1.369 
1.370 


0.983 
1.044 
1.096 


2.318 
2.459 
2.344 


0.0318 
0.0213 
0.0144 


67.65 
91.75 
115.7 


Navier- 
Stokes 


128 
256 
512 


0.0491 
0.0245 
0.0123 


2.502 
2.616 
2.545 


0.007 
0.00278 
0.00110 


1.482 
1.518 
1.551 


0.918 
1.268 
1.572 


3.549 
2.761 
2.376 


0.0247 
0.0114 
0.00539 


106.2 
150.4 
224.2 



Table 1: Parameters of the numerical simulations, number of collocation points N 3 ; grid 
spacing dx; time of fully developed turbulence tg; viscosity v; root-mean-square velocity 
Vrms = \Z2/3Eki n ; mean energy dissipation rate e; integral scale L = (2/3Eki n ) 3 / 2 /e; dissi- 
pation length scale rj = (i^/e) 1 / 4 ; Taylor- Reynolds number = \/ 15v rma L / v; all taken at 
the time of maximum enstrophy t = tg. 




Figure 1: Left: Comparison of the kinetic energy for Navier-Stokes, Burgers and the proposed 
model equation. The steep discontinuities of a Burgers flow explain the fast energy dissipation. 
Right: Evolution of free- falling turbulence for Burgers equation, the proposed model equation 
and Navier-Stokes equation. The turbulent flow is fully developed at t = tg, when the 
enstrophy € reaches its maximum. 



of Burgers equation but faster than the formation of vortex tubes for Navier- 
Stokes equation. Precise values for the timescales are stated in Table [T] 

Figure [2] depicts a volume render of the fully developed turbulence, the snap- 
shot in each case taken at t = tg. As expected, the Burgers flow (left) is domi- 
nated by shocks and the Navier-Stokes flow (right) consists of vortex filaments. 
For the proposed model equation (middle), the most dominant structures are 
two-dimensional folded vortex-sheets. 

A phcnomenological description, which takes into account the most dissi- 
pative structures of the turbulent flow, is the model of She and Leveque [H 
connected to log-Poisson statistics of the local energy dissipation Q. They 
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Figure 2: Dissipative structures in fully developed turbulence for different hydrodynamical 
models. Left: 2-dimensional shock fronts in Burgers turbulence. Middle: Folded vortex 
sheets in the new model equation. Right: Vortex filaments in Navier-Stokes turbulence. 



state that the scaling exponent £ p of the p-th structure function behaves like 



(5) 



where Co is the co-dimension of the most dissipative structures in the evolved 
flow and k is the time-scaling exponent. This formula will be referred to as the 
She-Leveque model even though in Q it is applied exclusively to the Navier- 
Stokes equation (where Co = 2, k = 2/3). 

The typical time U for the evolution of discontinuities for turbulent Burgers 
flows scales linearly with I, which accounts for k = 1. As the shocks traveling 
through the domain are two-dimensional we furthermore obtain Co = 1. Insert- 
ing this into equation © leads to ( p — 1. Since parts of the velocity field are 
continuous, for p < 1 the smoother regions of the velocity field are pronounced. 
This would equal a scaling exponent of Cp = P for these orders. Since this result 
is smaller than £ p = 1, it is dominant for p < 1. Therefore, we get 



p iorp < 1 
1 forp > 1 



(6) 



for Burgers equation. 



p 


1 


2 


4 


5 


6 


7 


Otp 


0.70 


0.95 


1.00 


1.00 


0.98 


1.00 


Cp 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 



Table 2: Scaling exponents ( p for Burgers equation. The a p measured with ESS are compared 
to the f p predicted by the She-Leveque model. 



Tabic [2] shows the measurements of a direct numerical simulation with 512 3 
grid points. Here, a p is the data obtained via ESS and £ p is the prediction of 
equation ([6]). The visualization of this result is shown in figure [3] on the left. 
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Figure 3: Scaling exponents for fully developed turbulence, comparison between theory and 
numerical simulation. Left: Burgers equation. Right: Presented model equation. 



Especially for high order of p (p > 3) the scaling exponent agrees with the 
prediction, yet smears out for smaller p (see also [J] ) . 

Similar to Burgers equation the introduced model equation allows for shocks 
to develop, since no incompressibility condition is stated. Nevertheless, energy 
decays solely via the dissipative term. Because of this, the hypotheses of Kol- 
mogorov may be adapted to Euler-Burgers equation. The Richardson cascade 
as well as the properties of the energy spectrum and the scaling behavior of 
structure functions should agree with the conclusions of Kolmogorov and the 
She-Leveque model. On the other hand, the structures that evolve seem to be 
significantly different from the vortex filaments known from Navier-Stokes. 

Figure [2] (middle) shows the most dissipative structures of fully developed 
turbulent flow. The norm of the vorticity visualizes two-dimensional folded 
vortex sheets as the structures that correspond to the vortex tubes of Navier- 
Stokes. This suggest Co = 1. The time-scaling exponent ti for the introduced 
equation is estimated as k = 2/3, with the same reasoning as for Navier-Stokes. 
Thus, 




is the prediction for the scaling exponents proposed by the She-Leveque model. 



p 


1 


2 


4 


5 


6 


7 


a p (Pu) 


0.4 


0.73 


1.22 


1.39 


1.55 


1.69 


(p 


0.418 


0.751 


1.213 


1.395 


1.556 


1.701 



Table 3: Scaling exponents f p for the new model equation. The a p for Pu measured with 
ESS are compared to the f p predicted by the She-Levque model. 

Table [3] features the results measured via ESS for the scaling exponents of 
Pu. As can be seen in figure [3] (right) the numerical data agrees very well with 
the prediction of the She-Leveque model. The solenoidal field is consistent with 
theory from low orders of p up to the highest order that was measured. 
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4. Global solutions 



In this section, we show global regularity for equation (01 for suitable initial 
data without any size restrictions. For this, we proof the remarkable property, 
that this equation possesses an infinite number of conserved quantities in the 
inviscid limit. Especially, the finitencss of the £ 4 -norm of the velocity field 
coupled with standard estimates for the ii^-norm enables us to show global 
regularity. 

The problem of whether the three-dimensional incompressible Navier-Stokes 
equations can develop a finite time singularity from smooth initial conditions 
or if it has global solutions remains unresolved (see 043 and the references 
therein). The answer to this important question is recognized as one of the 
Millennium prize problems [H, Q . 

Despite the complexity of the topic, a lot of progress has been made on this 
field in the past. For the two-dimensional case, global- in-time existence of unique 
weak and strong solutions is well-known (see |5, 6]). In three dimensions weak 
solutions are known to exist globally in time. For strong solutions, existence 
and uniqueness is known for a short interval of time which depends continuo usly 
on the initial data [l(J. Many results published in the past, starting with 
provide criteria for the global regularity of solutions via conditions applied to the 



velocity field [12|, [13( or components thereof [14| , the vorticity [15( , its direction 



iponents 

[161 or to thc pressure field [II, III]. 

The theory for the compressible Navier-Stokes equation is less well devel- 
oped, and we will not attempt a summary here. The multi-dimensional Burgers 
equation [l9| can be regarded as a crude simplification of this model. Global 
existence and uniqueness of strong solutions can be established in two and three- 
dimensions for suitably small initial conditions, much as with the Navier-Stokes 
system. Irrotational flows do possess global solutions for large data in arbitrary 



dimension, thanks to the Colc-Hopf transformation [2CJ, |2l|. However, there is 
no multidimensional weak theory because of the absence of a mechanism for 
energy dissipation, unlike Navier-Stokes. 

The situation for this new modified Navier- Stokes like equation is rather 
different. In this section, the existence of global solutions is proven for the model 
equation ((4]) in a domain ft which shall either be M 3 or a periodic cube in M 3 . 

Theorem 1. Let u G H 1 ^). Let f G Lf oc (R+,L 2 (Q)) n Lj oc (R+, L 4 (Q)). 
Then the initial value problem for the model equation has a unique global 
solution 

ueC(R + ,H 1 (il))nLl c (R+,H 2 (n)) . 

The aim is to show that the solution remains a priori bounded in L°°([0, T), i7 1 (fi))n 
L 2 ([0, T), H 2 (£l)) for any T > 0, which implies its existence and uniqueness with 
standard arguments comparable to e.g. fa. |22|. Throughout the argument, we 
denote the Euclidean norm of the vector u = J2i u i e i by u = (Si" 2 ) 1 ^- We 
first prove the following lemma: 

Lemma 1. Let Uo, f, O be defined as above. Then the quantity ||u(t)||j,p re- 
mains finite for 2 < p G R. 
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Proof. Taking the Euclidean inner product of with u yields the identity 

i (^-u 2 + Pu • Vii 2 ) = V -l\u 2 - ^Vu| 2 + f • u . (8) 



2 V<9i 

Integrate ((8]) over f2, use the fact that Pu is divergence free, and then apply the 
Cauchy-Schwarz inequality to obtain 

^l|u||! 2 +H|Vu|| 2 L2 <||f|| L2 ||u|| L2 . (9) 

Defining x(t) = \ (j|u(f)|| 2 2 + f Q ^||Vu(.s)|| 2 2 ds) , we have that 

x'(t)<\\f(t)\\L*(2x(t)) 1/2 ■ 
Upon integration, this gives the inequality 

\\u(t)\\ 2 L2 + j\\\Vu(s)\\ 2 L2 d s < (\\u \\ L , + J" \\f( S )\\ L ,d8\ . (10) 

With this estimate the lemma is shown for the case p = 2. 
Let 2 < n G R and multiply the identity © by u 2(,l_1) : 

1 fd 



± ^ + Pu • W«J = | (£-{u 2 ^)JL-u 2 ) - |Vu«| 2 ; 

-^("-^iVup + u 2 ^- 1 ^ u. 
Integrate this over il and apply Holder's inequality: 

|u||i'L + J |W| 2 + uu 2 ^)\Vn\ 2 ) d x < ||f|k-||u|| 2 L r 



2ndt' 



If we let 

y{t) = ^||u(t)|| 2 'L +f J |V W ™(.s)| 2 + ^ 2 (™- 1 )( S )|Vu( S )| 2 ) dxds , 

then we obtain 

y'(t)<||f(t)|| L2 „(2n 2/ (t)) 2 ^ i . 
This leads to the estimate 

||u(t)|| 2 'L + 2n £ J \Vu n (s)\ 2 + vu^-V^lVuis)^ dxds (11) 

< (||uo|| L -+/ 4 ||f( S )|| L ^) 2n . 

□ 



Remark: This key argument fails for the case of the Navier-Stokes equation. 
At the same time, this estimate establishes an infinite number of conserved 
quantity in the unforced inviscid case. 



Proof of Theorem [0 Take the L 2 -inner product of ((4]) with Au and integrate 
by parts: 



--I 

2 at 1 



Vu||| 2 +is\\Au\\ 2 L2 = / (Pu- Vu) • Audx+ / f-Audx 
Jn Jn 



(i) (it) 

The forcing term (ii) has the bound 



f f.Aurfx<||f|| L2 ||Au|| i2 < J||Auf L2 + i 
Jn 

The nonlinear term (i) is estimated as follows: 
^(Pu-Vu).Audx = _J^_ Ui ^_( (Pu)jg |_ Ui ) dx 

= -/^i(( Pu )iaf-^ + ^( Pu ^'4^) dx 



= 'I (KP^al-IVul 2 + ^«4(^(Pu) jM ,)) dx 

/ 7T~ 7n~ " ' 7U7 ' "< 
< !|V 2 u|| L2 ||VPu|| L 4||u|| L 4 . 

The second norm above is handled by interpolation. We first note that 

||VPu|| L , <||VPu||^ 4 ||VPu||^ 2 4 . 

Now when f2 = IR 3 , the Sobolev embedding theorem gives 

||VPu|| L6 <C||V 2 Pu|| L2 . (12) 

When n is a periodic domain, the norm on the right must be replaced by 
||VPu||jyi. However, since VPu has zero mean, this is bounded again by 
C||V 2 Pu|| £2 , by the Poincare inequality. Therefore, (|T2"j) holds in both cases. 
Using the facts that the operator P commutes with derivatives and that it is a 
projection in L 2 , we have that 

||VPu|| L2 < ||Vu|| L2 and ||V 2 Pu|| L2 < ||V 2 u|| i2 . 

Next, we use integration by parts to obtain the simple ellipticity identity 

_ r _d__d_ u ._d__d_ u . ^ x 

JQ dxj dxj % dxk dxk 1 

= l|Au|| 2 2 . 
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Combining these observations with Young's inequality, we conclude that the 
nonlinear term (i) is bounded by 




Altogether, we get the inequality 




Using Gronwall's inequality, we find that 

||Vu||i a ||Au( S )||i 2 < ||Vuo||| 2 exp^/ *||u(s)||| 4 d S (14) 

Jo 

+ £/o(^£/. t |M<0ll!< da) \\{(s)\\l 2 ds. 
Combining (fTU)) and Lemma Q] with p = 4, and (HU), we see that the quantity 

l|u(t)||iri+ / K||Vu( S )||| 2 + ||Au( S )||| 2 ]d S 
Jo 

remains finite. However, by CE3]) and the fact that Lf£ c (R + , i 2 (0)) C Lf oc (R+, L 2 (0)) 
we have that ^ 

\Ht)\\m+ I H|u(s)lllWi) ds 
Jo 

also remains finite. □ 

5. Final remarks 

In this paper, a modified Navier-Stokes equation is presented. Its dynamics 
and turbulent behavior are studied in terms of the scaling properties of its struc- 
ture functions. The most dissipativc structures arc identified as vortex sheets of 
co-dimension 1, which allows us to compare the numerically measured scaling 
exponents to a modified phcnomenologically based She-Lcveque approach. 

Furthermore, we prove the existence of global solutions for this equation. A 
remarkable consequence of this proof is the existence of an infinite number of 
conserved quantities ||u||^ P in the ideal (non- viscous) case without forcing. This 
property is not only responsible for the existence of global solutions but should 
show up in the statistics of intermittent turbulent fluctuations. Work in this 
direction is in progress. 
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